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We use holography to study Conformal Phase Transitions, which are believed to be realized 
in four dimensional QCD and play an important role in walking technicolor models of 
electroweak symmetry breaking. At strong coupling they can be modeled by the non- 
linear dynamics of a tachyonic scalar field with mass close to the Breitenlohner-Freedman 
bound in anti de Sitter spacetime. Taking the action for this field to have a Tachyon-Dirac- 
Born-Infeld form gives rise to models that resemble hard and soft wall AdS/QCD, with a 
dynamically generated wall. For hard wall models, the highly excited spectrum has the KK 
form nin ~ n; in the soft wall case we exhibit potentials with m„, ~ n", < a < 1/2. We 
investigate the finite temperature phase structure and find first or second order symmetry 
restoration transitions, depending on the behavior of the potential near the origin of field 
space. 



1. Introduction 

Consider SU{N) gauge tlieory coupled to F massless Dirac fermions in the funda- 
mental representation of the gauge group.0 For F > llA^/2, this theory is free at long 
distances, and its low energy dynamics can be studied using the techniques of perturbative 
field theory. For F < llN/2, the theory becomes interacting in the infrared where it is 
said to be in a non-abelian Coulomb phase. The crossover between the free ultraviolet 
behavior and the interacting infrared conformal field theory occurs at a scale Aqcd- As 
F decreases, the infrared coupling increases; eventually, perturbation theory fails at scales 
of order Aqcd and below. 

On the other hand, when F N the theory is believed to confine and break the chiral 
symmetry SU (F) l x SU (-F)_b — )■ SU (-F)diag- The 't Hooft large analysis of perturbation 
theory provides a nice picture of what happens in this regime. To leading order in F/N, 
the only fields that run in loops are the adjoint degrees of freedom (the gauge fields) . The 
spectrum includes — 1 massless Nambu-Goldstone bosons ("pions"), and a discrete 
spectrum of massive glueballs and mesons with masses of order Aqcd- 
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Fig. 1: Schematic phase diagram of QCD as a function of tlie number of flavors 
F. Aqcd is tlie scale above which the theory becomes free; /i is the meson mass 
scale. 

It is natural to ask how the above two pictures are connected as one varies the number 
of fiavors F. Since the realization of chiral symmetry is different in the two regimes, there 
must be a phase transition at a finite value of F/N. This transition is believed to be 
continuous; the order parameter, which can be taken to be the typical meson mass /i, is 

^ It is convenient for our purposes to study this theory in the limit F,N ^ 1, so that the 
parameter F/N can be treated as continuous. Much of what we say can be extended to finite 
F,N. 
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expected to behave as a function of F in the way depicted in figure 1. According to this 
picture, the non-abehan Coulomb phase extends down to a critical number of flavors, F^, 
which is not known precisely but is believed to be Fc — 4:N. When F decreases past 
Fc, n becomes non-zero, but near the transition it is much smaller than the scale Aqcd 
above which the theory becomes free. As F decreases further, n increases and eventually 
approaches Aqqd- 

The transition at F = Fc is not well understood. A widely discussed scenario (see 
e.g. and references therein) is that it is driven by the gauge invariant fermion bilinear 
operator tl^t/;. The UV scaling dimension of this operator is three, but in the non-abelian 
Coulomb phase its IR scaling dimension is lowered by gauge interactions. It is believed that 
as -F — 7- Fc, the IR dimension approaches two, and when F < Fc it becomes complex. This 
behavior signals an instability of the Coulomb phase due to a non-vanishing /3-function 
for the "double trace" operator (Tr-^-i/')^, which leads to the appearance of the small mass 
gap fi in figure 1. Near Fc, the gap scales like fx ~ AQCDQ^pi—a/VF — Fc), which is 
usually referred to as Miransky scaling. The corresponding phase transition is known as a 
Conformal Phase Transition (CPT) f^. 

Understanding the behavior of the theory in the vicinity of the phase transition at 
F = Fc is an interesting open problem, which may also be important for applications. In 
particular, a version of technicolor known as "walking technicolor" relies on it to gener- 
ate electroweak symmetry breaking [^, and similar dynamics may have applications to 
condensed matter systems such as graphene (see e.g. for a recent discussion). 

In thinking about the physics near the transition, it is useful to divide the RG flow of 
QCD into two stages. The flrst involves the flow from an asymptotically free fleld theory at 
energies well above Aqcd to a strongly coupled CFT well below that scale. The only part 
of this flow which is of interest for studying the transition is its endpoint, the non-trivial 
infrared CFT which serves as the arena for the transition. Therefore, for our purposes we 
can treat Aqcd as a UV cutoff, and focus on the dynamics well below it. For F > Fc, the 
infrared CFT is stable under RG, while for F < Fc it dynamically generates the scale /x. 
The problem one is faced with is to understand the mechanism of mass generation. 

As stressed in , the behavior exhibited by this example occurs much more generally, 
whenever two flxed points of the RG approach each other, merge and move off into the 
complex plane. In general (at large A^), the transition is driven by an operator O whose 
dimension approaches d/2 near the critical point. In QCD, d = 4 and O = Tr'0'0; examples 
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with other values of (i, which involve systems of fermions localized on defects interacting 
with gauge fields in the bulk were recently discussed in 

The fact that the transition is continuous makes it natural to expect that it can be 
described by studying the dynamics of the order parameter O. A conventional Landau- 
Ginzburg description is inappropriate here, since the transition is essentially an infinite 
order one (as reflected e.g. in the Miransky scaling). A further difficulty is that the 
dynamics is typically strongly coupled in the vicinity of the transition. It is natural to 
ask whether holography can be used to shed light on this problem. In this paper, we will 
take some steps towards this goal. Other recent discussions of holographic walking include 

The starting point of our discussion is the assumption that the d dimensional CFT that 
serves as the arena in which the phase transition occurs has an AdSd+i dual. The order 
parameter O corresponds in the bulk description to a scalar field T whose mass approaches 
the Breitenlohner-Freedman (BF) bound. Thus, to describe the physics associated with 
the scale fx in figure 1, we need to study the dynamics of T in AdSd+i- 

Using the standard 't Hooft map [^j, the field T can often be thought of as an open 
string tachyon.il The dynamics of such tachyons has been extensively studied in the past 
and a lot is known about it (see e.g. []T^ for a review). One of the interesting outcomes of 
these studies is the realization that the dynamics of open string tachyons is well described 
qualitatively, and in some cases quantitatively, by an effective action called the Tachyon- 
Dirac-Born-Infeld (TDBI) action [|T7| - |T9| . Thus, it is natural to model the phase transition 
using such an action in AdSd+i- 

The TDBI action depends on the choice of a potential V(T) for the tachyon. The 
quadratic term in the potential gives the mass of the tachyon, which we will take to be close 
to the BF bound. As we will see, the low energy dynamics depends on the full potential, 
and one of our main goals will be to explore this dependence. This can be viewed as a 
bottom-up approach, analogous to the one taken in AdS/QCD. We will see that when the 
mass of the tachyon is slightly above the BF bound, the theory is in a conformal phase, but 
when the mass is slightly below the bound, it breaks conformal symmetry and generates a 
mass gap. This gap can be made parametrically small relative to the UV cutoff A, which 
plays the role of Aqcd- Depending on the choice of potential, we find different spectra of 
mesons, including some that agree with what is expected in QCD in the 't Hooft limit. 



^ This is certainly true in probe brane systems, such as those studied in and formally is 
also true in QCD, although for F ^ N the distinction between open and closed strings is blurred. 
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We will also explore the thermodynamics of TDBI models and find that they exhibit 
first or second order phase transitions at a finite temperature, depending on the small T 
behavior of the potential V{T). This difi'ers from most other holographic models, which 
exhibit strongly first order phase transitions. We will discuss the origin of the differences 
between our model and other models of holographic QCD in section 7. 

The plan of the paper is as follows. In section 2 we revisit a system that we studied in 
, = 4 SYM coupled to fermions in the fundamental representation of the gauge group 
U{N) localized on a (i — 1 dimensional defect. There, we showed that this system exhibits 
a conformal phase transition, and studied it using holography. In section 2 we show that 
the discussion of can be rephrased in terms of a TDBI action for a real tachyon field T, 
and we compute the potential V{T). This provides support for the idea that one can use 
the TDBI action to describe CPTs such as the one expected in QCD. 

In section 3 we discuss the TDBI approach to general CPTs. We describe the vacuum 
solution for different tachyon potentials, and in particular show that the TDBI model gives 
an interesting dynamical realization of soft and hard wall models of AdS/QCD. 



We also derive the Lagrangian for small fiuctuations around the vacuum, which correspond 
to scalar and vector mesons. 

In sections 4 and 5 we discuss the hard and soft wall models, respectively. We show 
that hard wall models always have KK-type spectra of excitations mn ~ n at large excita- 
tion number n, in agreement with the analysis of which involved a particular hard wall 
model. On the other hand, soft wall models can give other types of asymptotic spectra. 
We discuss explicitly a class of potentials which gives the asymptotic spectrum ~ n" 
with arbitrary < a < 1/2. We also discuss the question of restoration of chiral symmetry 
in the asymptotic spectrum of (axial) vectors and (pseudo) scalars. 

In section 6 we study the thermodynamics of TDBI models. We show that they 
undergo a first or second order phase transition at a finite temperature, depending on the 
form of the potential near the origin of field space. The structure we find is rather different 
than in other holographic models of confinement. In section 7 we discuss our results and 
comment on possible extensions. Some technical details of the meson spectrum in the 
axial sector, and a numerical check of the finite temperature analysis, are relegated to the 
appendices. 



2. Defect fermions coupled to = 4 SYM 

In [0 we showed that the system of A/" = 4 SYM coupled to fermions localized on a 
d — 1 dimensional defect exhibits a transition similar to that expected in QCD. While in 
QCD this transition occurs as a function of the number of colors and flavors, in the system 
studied in the transition happens for any number of flavors, and is driven by the 't Hooft 
coupling A of A/" = 4 SYM. As the coupling increases, the scaling dimension of the fermion 
bilinear -i/^-i/; decreases. When the coupling reaches a critical value Ac, A{i(^i(^) approaches 
(i/2, and the system undergoes a transition from a conformal phase to one in which the 
fermions are massive. This is an example of what we referred to in the introduction as a 
conformal phase transition (CPT). 

The defect theory can be embedded in string theory as the low energy theory on 
N iD3-branes and one or more Dp-hranes, which intersect the DS-branes on a d dimen- 
sional spacetime and are extended in n additional directions transverse to the threebranes. 
Strings stretched between the D3 and Dp-hranes give rise to fermions in the fundamental 
representation of U{N). At weak coupling, the dynamics is described by a Lagrangian 
which couples the fermions to the DS-brane gauge fleld and 6 — n = d — p + 5 adjoint 
scalars. At strong coupling, one can replace the Z)3-branes by their near-horizon geometry 
AdS^ X S^, and study the dynamics of the L'p-branes as probes in this geometry. 

In we pointed out that one can study the defect system as a function of the 
parameters d, n, formally viewing them as continuous variables, in the spirit of the e 
expansion. By varying these parameters, one can change the critical coupling Ac- For flxed 
n and d slightly above two, Ac ^ 1 and the transition can be studied using perturbative 
gauge theory. For flxed d > 2 and n slightly above 1 -|- (i^/4, Ac is large, and the transition 
can be studied using holography. At strong coupling, the transition is between a phase in 
which the probe brane wraps an AdSd+i x S^~^ inside AdS^ x S^, and one in which the 
probe brane is deformed in the IR. In we used the gravitational description to study 
aspects of the physics associated with the transition. 

The purpose of this paper is to generalize the discussion of [0 to a larger class of 
systems. The eventual hope is to construct a holographic description of QCD in the 
vicinity of the transition from non-abelian Coulomb to conflning infrared behavior. As a 
step in this direction, we will present a class of bottom up models whose properties are in 
qualitative agreement with expectations. As we will see, our construction makes contact 



with both the hard and soft wall models of AdS / QCD [pO| , |2T[] , although the way they arise 
here is rather different than in the literature. 



5 



We start with a brief review of some aspects of . The near-horizon geometry of the 
Z)3-branes, AdS^ x 5"^, is described by the metric 




where x^, /i = 0,1,2,3, are directions along the L>3-branes, and r and O5 are spherical 
coordinates on the transverse IR^. L is the radius of curvature of AdS^ and the size of the 



Into this background, we place a Dp-hrane that wraps d of the four x^ and an IR" 
subspace of the transverse IR^. Using spherical coordinates (p, fln-i) on IR"^, we can write 
the metric ( ^.11 ) as 



The Dp-hrane wraps an AdSd+i x 5""^"^ parametrized by x", a = 0, 1, 2, ■ ■ ■ , (i — 1, as well 
as p and 0„_i. The coordinates c^^, i = 1, ■ ■ ■ , 6 — n, parametrize directions transverse to 
the probe brane. The radial coordinate r is given by 



Excitations of the probe brane correspond to fluctuations of the scalar fields (pi, the world- 
volume gauge field, and fermions. Their dynamics, which is governed by the DBI action, 
was studied in 0. 

It is useful to note that the action of the probe brane must be SO{d, 2) invariant. 
The reason is that the action can be thought of as describing small fluctuations around the 
AdSd+i y<S^~^ configuration (pi = 0. Although this configuration is not always dynamically 
stable, this does not matter for the symmetry structure. SO{d,2) is the subgroup of the 
50(4,2) isometry of the metric ( |2.1| ), which acts on the coordinates (r, x"). 

In the DBI action derived from the metric ( |2.2|) , which was discussed in 0, the above 
SO{d, 2) symmetry is hidden for the following reason. While the symmetry acts naturally 
on the radial coordinate r, we have split it here as in ( p.3|) , with p a worldvolume coordinate 
and (pi fields on the brane. It is thus clear how one needs to modify the analysis of 
to make the symmetry manifest .i The idea is to take the probe Dp-hrane to wrap the d 

^ We thank E. Martinec for a useful discussion of this issue. 




(2.2) 




(2.3) 
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dimensions labeled by x°', the radial direction r, and an S'^~^ of maximal radius in the 
S^. The 6 — n directions on transverse to the S'^~^ give rise to scalar fields on the 
worldvolume. These scalar fields are by construction compact, since they correspond to 
angular variables on the sphere. Their action is manifestly invariant under SO{d,2). 

To make this more explicit, consider an excitation of one of the scalar fields 9 = 9{x^), 
where (x^^) = (a;",r) parametrize the AdSd+i- The metric on the can be written as 



dn'i = L'^de^ + {L cos 9Ydni_^ + • • • (2.4) 

where the ellipsis denotes terms associated with other directions on the sphere, which we 
do not excite. The induced metric on the Dp-brane now takes the form 

ds"^ = [qmn + L^dMOdNe) dx^dx^ + (L cos efdnl_^ (2.5) 



where Qmn is the AdSd+i metric ((pTll) restricted to x^). The DBI action that follows 
from (p75| ) is@ 

S = - j d'^+^x^/^{cos ey-^ + L^g^^dMOdNO. (2.6) 

This action is manifestly covariant in AdSd+i and thus is SO{d,2) symmetric. From the 
point of view of the parametrization employed in 0, the excitation 9{x^) corresponds to 
one of the 6 — n scalar fields, say (f) = (f)i. The map between the coordinates (p, cf)) of 
and (r, 6) here is 

p = rcos9; (f) = rsm6. (2-7) 
The action (|2.6|) has the TDBI form that was studied in the context of open string 



tachyon condensation [16-19,22-25 



S^- j d'^+^xViT)y/^ ^ - J d'^+^Xy/^ViT)^! + g^^duTdNT, (2.8) 

where 

Gmn = 9mn + duTd^T (2.9) 
is the induced metric, G = det Gmn , T = L9 is the "tachyon" field, and 

rp\ n—1 



y(T)= (cos-) (2.10) 



We omitted an overall multiplicative factor. 
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is the tachyon potential. Note that the tachyon field is bounded, \T\ < Tjji = Ln/2. The 
potential ( |2.10|) vanishes at the edges of field space, a fact that will play an important role 
in our discussion. The mass of the tachyon, obtained by expanding (|2.10| ) to quadratic 
order and writing V{T) = 1 + ^m^T^ + ■ ■ ■, is 

m' = -^. (2.11) 

As a check, according to the AdS/CFT dictionary, the asymptotic behavior of a Klein- 
Gordon field of mass m is r~^, where A(A — d) = 7v?L^. In we found the large r (or 
small T) behavior T ~ j.-d/2±iy/^^ with k = n — 1 — d'^ /4, which agrees with ( 2.11 ). 

To summarize, we see that the discussion of [0 can be phrased in terms of the TDBI 
action ( |2.8| ). The microscopic model fixes the tachyon potential ( p.lO| ). The potential 
varies between one at T = 0, which corresponds to the Dp brane located at = in the 
parametrization (|2.2|), and zero at T = Tir. 

By mapping the analysis of to the tachyon DEI language, we can deduce several 
properties of the Lagrangian ( |2.8|) . The vacuum of the theory depends on whether the 
mass of the tachyon, ( |2.11| ), is above or below the BF bound. For 

m2>m|^ = -^, (2.12) 

the vacuum corresponds to T = and preserves conformal symmetry. For smaller m^, the 
vacuum is described by a non-trivial T{r) which can be obtained by minimizing the energy 
(P78|). To do this, one has to introduce a UV cutoff. In 0] the cutoff was introduced by 
picking a maximal value of p < A and setting (p{p = A) = 0. In the TDBI coordinates, 
this corresponds to setting r < A and imposing the UV boundary condition T(r = A) = 0. 

Since the equation of motion for T(r) that follows from ( |2.8| ) is second order, we need 
one more boundary condition to fix the solution. In we required that the brane be 
smooth at small p, which led to dp(p{p = 0) = 0. With these boundary conditions, we 
found that the lowest energy state of the brane has the shape depicted in figure 2. This 
configuration breaks conformal symmetry and dynamically generates the mass scale /i. /i 
can be computed in terms of the UV cutoff and the distance from the transition k; it 
satisfies Miransky scaling. One can take the double scaling limit A — t- oo, k — )■ with p 
held fixed, which was referred to in as the BKT limit (using the terminology of 0). 
This limit focuses on the physics of the conformal phase transition while decoupling the 
dynamics associated with the UV cutoff. 



8 



A 



P 



Fig. 2: The shape of the brane parametrized in terms of (p, t/)) and {r,9). 

T . 

T ~ 

IR 



|j, A r 

Fig. 3: The shape of the brane in figure 2 parametrized in terms of T(r). 

In terms of the TDBI variables T, r, the shape of the brane is depicted in figure 3. 
Note that: 

(1) In the trivial vacuum T(r) = which preserves conformal symmetry, r ranges from 
zero to infinity. In contrast, in the vacuum with broken symmetry the radial coordinate 
is bounded from below {r > n). This is reminiscent of what happens in holographic 
models of confinement, but there the excision of a finite region in r is usually put in 
by hand, as in the hard wall model of AdS/QCD [|0], or arises due to closed string 
(gravitational) dynamics [p6| , |27[| . Here, the effect arises due to the dynamics of IR 
unstable open string modes.i 

(2) In the (p, (j)) parametrization used in 0, the fact that the dynamically generated mass 
fi could be made arbitrarily small by tuning k meant that the field (^(p) describing the 
vacuum with broken symmetry could be made arbitrarily small everywhere (see figure 
2). In the (r, T) parametrization the picture is somewhat different. No matter how 
small /i is, the field T(r) always interpolates between T = at r = A and T = Tm 
(recall that Tm = Ln/2 is the zero of V{T)) at r = p. The scale p controls the size 
of the excised region in r space, rather than the range of values taken by T(r). Of 
course, as p decreases, the profile T(r) remains small in a larger range of r. 

^ See also comment (4) below. 
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(3) The IR boundary condition (f)'{p = 0) =0 translates in the (r, T) parametrization to 
the condition that r(T) is stationary at Tm, 

^{T = Tin)=0, (2.13) 

or equivalently, drT diverges as T — )■ Tm. r{T, x"), thought of as a field living on the 
interval < T < Tm, satisfies Dirichlet boundary conditions at T = and Neumann 
ones at T = Tjfi. Thus, the scale p = riTm), which governs the size of the region in 
r that is excised from the space, is dynamical. 

(4) The preceding discussion is reminiscent of what happens in the holographic models 
of spontaneous symmetry breaking studied in [^-^| following [31 1. Indeed, if we 



rotate figure 3 by ninety degrees, so that the T and r axes become horizontal and 
vertical respectively, we arrive at a similar picture to those papers. The symmetric 
configuration T(r) = can be thought of as describing a brane extending along the 
r axis from r = A to r = 0, then going up the T axis, from T = to T = T//^. 
For rm? < rn\p this configuration is unstable and is dynamically deformed to that of 
figure 3, which looks like half of the U-shape describing the symmetry breaking phase 
28| - p0[] . The boundary condition ( |2.13| ) is very natural from this point of view. It 



m 



ensures that one can continue the configuration of figure 3 past T/i^ and connect it to 
its mirror image, obtained by refiecting the configuration of figure 3 about the dashed 
line. This gives a smooth U-shape, very similar to those found in the above papers. Of 
course, there are important differences between the systems. In particular, comparing 
the action (|2.8[) to the one given by eq. (2.11) in PO], we see that the analog of V(T) is 



constant there, while the function of r in front of the square root, which here is r'^"^, 
is in general different there. This accounts for the different physics of these systems. 
An interesting question is whether the analysis of this section can be extended to more 
general CFTs which generate a mass scale due to the dynamics of operators with dimension 
close to ci/2. In holography, such operators correspond to scalars with mass close to the 
BP bound, and if they come from the open string sector, it is natural to try to describe 
them by a tachyon DBI action with some V{T). In the next section we will develop such 
a description, and study the resulting models for different V{T). 
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3. TDBI description of dynamical symmetry breaking 

Consider a (i-dimensional CFT which contains an operator with scaling dimen- 

sion close to d/2. If this CFT has a holographic description, the operator O is dual to a 
scalar field T(x^) in AdSd+i, whose mass is close to the BP bound (|2.12| ). Suppose that 
the mass depends on a tunable parameter. Then it may happen that as the parameter is 
varied, the mass approaches the BF bound and crosses it. In the CFT, this corresponds 
to the dimension of O approaching d/2 and going into the complex plane. While such 
behavior may appear strange at first sight, it is realized in the system studied in and 
is believed to be realized in (large A^) QCD. It was also discussed in the context of the 
AdS/CFT correspondence in ||3^ - |35|] and in the presence of background fields that break 
Lorentz symmetry in [p6| - ^ . 



In the systems studied in []32| - |35[| the bulk field T dual to (9 is a closed string tachyon, 
while in the others it comes from the open string sector. As mentioned in the introduction, 
some aspects of the dynamics of open string tachyons are well described qualitatively, and 
in some cases quantitatively, by the tachyon DBI action. It is thus natural to ask whether 
the TDBI action can describe the dynamics associated with more general CPTs than the 
one discussed in and the previous section. In this and the following sections, we will 
study such models. 

In the system discussed in the previous section, the "tachyon" was a geometric mode, 
and the "tachyon DBI" action was just the usual DBI action for the brane system. In the 
literature on open string tachyon dynamics, there have been proposals that one can more 
generally interpret the open string tachyon geometrically (see e.g. p3| - |25[] ). Although such 
an interpretation fits naturally with our results, we will not assume it here. 

The starting point of our discussion is the TDBI action ( p.8| ). This action describes 
a real tachyon field T on AdSd+i, whose dynamics depends on the choice of a potential 
V{T). We will assume that the field has a mass slightly below the BF bound, i.e. the 
potential V{T) has a local maximum at the origin: 

F(T) = 1 + ^m^T^-f... (3.1) 

where m? < m?^p = — (i^/4.i We will further assume that V{T) is a monotonically 
decreasing function of T that goes to zero at a particular value of T, T = Tjn. As we 



Here and below we take the AdS radius to be L = 1. 
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will see shortly, there is a qualitative difference between the cases of finite and infinite 
TiR (see figure 4). In the system of Tir was finite, while for the potentials that arise 
in open string tachyon condensation it is typically infinite P3| , ^ . We will consider both 
cases below. 



V 










T 




Fig. 4: The tachyon potential, for finite (a) and infinite (b) Tir. 



The fact that the potential vanishes somewhere in field space is important. In the 
study of open string tachyon condensation this value of the tachyon corresponds to the 
state in which the unstable D-brane (or brane-antibrane pair) disappears [jl^ . In the defect 
fermion system of it is the state where the probe D-brane is pushed to infinite (p, so that 
the fermions have an infinite current mass and decouple. Such a decoupling limit should 
exist for other systems described by ( |2.8| ) as well. 

One can think of the TDBI action as an analog of the Landau-Ginzburg action for 
an order parameter in the vicinity of a first or second order phase transition. Here, the 
order parameter T is a bulk field, which describes an infinite number of light states in 
the boundary theory. The potential V(T) labels different "universality classes" of such 
transitions. 

Before turning to a more detailed investigation of the action ( p.8| ), we comment on 
some generalizations. If the boundary CFT we are studying has a U{1) global symmetry 
(under which the operator O is not charged), generated by a conserved current J", we need 
to include in the bulk description a gauge field A^^ . This can be achieved by replacing 
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Gmn Gmn + Fmn in (|2.8| ). The modified action can be used to calculate correlation 
functions of the current J" and the spectrum of vector mesons, and to turn on a chemical 
potential for the corresponding conserved charge. 

In QCD, the tachyon field is complex, like the dual fermion bilinear This can 

be accomodated by using the TDBI action for the D — D tachyon TOJT^] , 



(3.2) 



where 



G^MN = 9MN + ^D(^mT^Dn)T + Fjil (3.3) 



and similarly for L ^ R. In QCD with Nf flavors, the tachyon transforms in the (A^/, Nf) 
representation of the SU{Nf)L x SU{Nf)ji global symmetry group. DBI-type actions for 
matrix fields are known to be inherently ambiguous. For the purpose of our discussion, it is 
enough to study the full non-perturbative action ( |3.2D for the trace of the Nf x Nf matrix T, 
which develops a non-trivial vacuum expectation value. The rest of the components of T, as 

(L R) 

well as the SU{Nf) l x SU (A^/)i? gauge fields A)^' which are dual to the global symmetry 
currents, can be treated perturbatively. In particular, to study the spectrum of mesons we 
only need to work to quadratic order in these fields, in which ( p.2|) is unambiguous. 

The vacuum of the model ( p.8|) is described by specifying a field configuration T = 
To(r) which minimizes the energy function obtained from (I 



£ = J drr'^-W{T).yi + r^T'{r)^. (3.4) 



The equation of motion following from ( |3.4| ) can be written as 

r'^-^drlnV _ d / r^+^T' \ 

or equivalently, 

-{1 + r^T'^)dT\nV + {1 + d)rT' + dr^T'^ + r^T" = 0. (3.6) 

We believe, but have not proven in general, that for any potential of the qualitative form 
in figure 4, the vacuum is trivial {i.e. To(r) = for all r) when the mass (|3.1|) is above the 
BF bound. We will mention some evidence supporting this claim below. We will also see 
that for m'^ below the BF bound the vacuum corresponds to a non-trivial solution, To(r), 
which takes the qualitative form depicted in figure 3. The boundary conditions satisfied 
by the solution are the same as those of section 2: T(r = A) = in the UV, and (|2.13|) in 
the IR. For models where Tm — )■ oo, there is no excised region and To(r) takes the form 
shown in figure 5. It can be thought of as a limit of the solution in figure 3, obtained by 
taking Tm — )■ oo and /U — i- 0. 
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Fig. 5: The vacuum configuration To(r) for infinite Tir. 

As is clear from figures 3 and 5, at large r the solution of (|3.6| ) is small, and we can 
replace the full equation of motion by the linearized one, 

r^T" + {1 + d)rT' - m^T = 0. (3.7) 

The general solution of this equation is 

To(r) = A(^)%in(^v^ln^ + (/)^ (3.8) 

where 

K = m\p — (3.9) 

is assumed to be positive, /i is a scale, and A, 4> are dimensionless constants whose definition 
depends on ii. As r decreases, To(r) increases, until eventually the linear approximation 
( P77| ) breaks down and one has to go back to the full equation of motion, n can be thought 
of as the scale at which this happens. The advantage of this definition is that it applies 
both to models with finite Tjn, where here is of the order of the scale /U defined in figure 
3, and with infinite Tjr, where there is no excised region. Of course, this definition is not 
precise, but it can be made precise e.g. by setting n to be the mass of the lowest lying 
meson. Other definitions differ from this one by numerical factors. The ambiguity in the 
definition of the dynamically generated scale is familiar in QFT, and it can be dealt with 
using the renormalization group. 

In the BKT limit A — )■ oo, k — )■ with held fixed, the large r behavior of T'o(r) 
becomes 

To(r)= (^)^ {c^Xrv-^+C^ (3.10) 
14 



where Ci, C2 are again dimensionless constants, related to A, cp as follows: 



Ci = Av^; C2 = A(j). (3.11) 

The presence of a non-trivial tachyon condensate To(r) leads to a deformation of the 
AdSd+i metric Qmn to 

Gmn = Qmn + dMTodNTo (3-12) 

and to a non-trivial dilaton $(r). Excitations described by the TDBI action live in this 
modified, or open string, background. As r — )■ 00, the open string metric G approaches 
the closed string (AdS) metric g and the dilaton approaches a constant, but for small r, 
the open and closed string backgrounds are in general quite different. In the rest of this 
section, we will derive the equations describing fluctuations of the TDBI action around the 
vacuum T'o(r). 

We start with fluctuations of the fleld T{x^), which give rise to u-mesons. We deflne 
the fluctuation y via T(x^) = To(r) +y{x^), plug this into ( p.8|) , and expand to quadratic 
order in y. The resulting action for the fluctuations is 

S2 = -^ J d'^xdrV^e-'^ {G^^^dMydNV + m^{r)y^) (3.13) 

where 

V{To) 



e 



2, , VTT^ 

m (r) =- 



V 



d'^V 1 ^^0 dV d ( r'^Tl^ 

[a — ij 



9T2 ^1 + r2T^2 ^ ^ dT + r^T^2 QT dr ^ ^\ + r'^T^ 

(3.14) 

and Gmn is the open string metric (|3.12| ). It is easy to check using (|3.1| ), ( p.lOj) , that 



at large r the dilaton $ approaches a constant, and the mass function m'^{r) goes to the 
tachyon mass m^. 

To study the spectrum of mesons, it is convenient to change coordinates from r to 
z(r), deflned by the requirement that the open string metric Gmn takes the form 

GMNdx^ dx^ = hapdx'^dx^ = r'^{z)dx"'dx^rjo,p. (3.15) 



Here (x") = (2, x"). Comparing (|3.12[ ) and ( |3.15| ) we see that the coordinates z and r are 
related as follows: 

^ = .^Upl. (3.16) 

dr 
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For large r, ( |3.16 ) reduces to z = 1/r, and ( |3.15| ) becomes a familiar parametrization of 



the metric on AdSd+i- For small r, there are significant deviations from AdS, which we 
will discuss below. 

The quadratic action for fluctuations ( |3.13| ) is given in these coordinates by 



= -\j d'^+^xr'^-^e-'^ (d^yd'^y + m^iryy^) 



(3.17) 



where on the second line r is viewed as a function of z, and indices are raised and lowered 
with the flat metric rjct/s. We next deflne the wavefunction ip = ye~^^'^ where 

B{z) = ^{z)-{d-l)\nr{z), (3.18) 

and write it as il){x"') = il;{z)exp{ikaX"'). The equation of motion for il;{z) that follows 
from ( |3.17| ) takes the Schrodinger form 

-ij" + Ve«{z)ij = M^ij (3.19) 

where the prime denotes differentiation w.r.t. z, Ai"^ = —kak"' is the meson mass, and the 
effective potential is 

V,s{z) = -^{B'f - ^B" + r{z)V{riz)). (3.20) 

The large z (IR) behavior of the potential Veff{z) depends on the details of the tachyon 
potential V{T). For small z, however, the potential is universal: To(r) goes to zero as 
r — 7- oo, the dilaton $(r) and mass function m(r) approach constants, and r = 1/z. Thus, 
in this regime one flnds B{z) = {d — l)lnz and 

Vesiz) = ( + = - f7 + ^H- (3-21) 



We are interested in k > (the massive phase), in the presence of a UV cutoff, which 
provides a lower bound z > zuv = 1/A on 2;. The boundary condition there is "ipi^zuv) = 0. 
If the potential 14fF was given by ( |3.21| ) for all z > zjjv, the spectrum of ( |3.19[ ) would 
contain normalizable states with M.^ < 0, i.e. tachyons. The mass squared of the lowest 
lying of these tachyons is of order — /x^, where /i ~ Aexp(— yr/y^). It is easy to understand 
this from the point of view of the TDBI system. The potential (|3.21| ) describes small 
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fluctuations about the conformally invariant state T(r) = 0. This state is unstable to 
condensation of T{r) towards the configuration To(r) described in figures 3 and 5. The 
tachyons of ( p.21| ) are the eigenmodes of this condensation process. 

When studying perturbations around the stable configuration T'o(r), the potential 
takes the form ( |3.21| ) for z ^ but is deformed at z ~ and larger. Since the 



wavefunctions corresponding to the tachyons of ( |3.19| ), ( |3.21| ) have finite support at z ~ 
the spectrum of Ves in the state corresponding to To(r) is deformed relative to that of 
the configuration T(r) = (even in the BKT limit). In particular, it cannot have tachyonic 
modes by construction, since it is the lowest energy state with the prescribed boundary 
conditions. 

Notice that for k < 0, the spectrum of (|3.19|) , ( |3.21| ) does not contain any bound 



states which would correspond to tachyons. This supports our earlier assertion that in 
the conformal phase, the state T(r) = is the ground state of the system - it shows that 
T(r) = is at least locally stable. The fact that k = lies on the boundary between 
two different behaviors of the bound state problem corresponding to the potential ( |3.21J ) 
is familiar from the QM analysis of this potential (see e.g. |]45[] ). 

To study the spectrum of vector mesons, we have to add a gauge field Am to the 
TDBI action. The action then becomes 

^ = -J d''+^xViT)^-det{GMN + Fmn) ^ So-^ j d''+'x^/^V{T)G^'^' G''''' FmnFm'N' 

(3.22) 

Thus, the gauge field propagates in a spacetime with metric Gmn ( |3.12D , and dilaton 

e""^ = V{To). (3.23) 

It is interesting that while the metric that the vectors and scalars experience is the same, 
the dilaton is not - it is given by (|3.14 ) for scalars and ( |3.23|) for vectors. Another difference 



is that the analog of m(r) in ( |3.13D vanishes for the gauge field. 

It is convenient to change coordinates as in ( p.l5|) , (|3.16| ), choose the gauge = 0, 
and write the remaining components of the gauge field as Aa = ^a(a^'')exp(i?/2)'0(z), where 
is the (transverse) polarization of a gauge field of momentum k"". The wavefunction 
i^iz) then satisfies a Schrodinger-type equation of the form ( |3.19[ ). 14ff is given by ( |3.2C1| ) 
with m(r) = and 

B{z) = -lnV(To) - {d-3)\nr{z). (3.24) 
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The small z (UV) behavior of the potential 14fj is again insensitive to the details of the 
tachyon potential; it is 14ff = [d — \){d — 3)/(42^). Thus we see that in the UV, the 
potentials for the a and p mesons both go like = —a/z^, but while takes the critical 
value 1/4, ctp is strictly below this value. This may be related to the fact that while the 
(T- meson spectrum was found in to contain a light mode which can be thought of as 
a pseudo-Nambu-Goldstone meson of broken conformal symmetry, the p-mesons were all 
heavier. 

One also needs to check whether fluctuations of the gauge field Am give rise to scalars. 
In the gauge = 0, the scalar modes may be parametrized as Aaiz.x"-) = 4>{z)daT^{x^). 
It is easy to check that the equations of motion of ( |3.22|) do not allow such solutions with 
non-zero mass. The case of zero mass was analyzed in [^, where it was shown that a 
massless meson appears if the integral 

^^ = j dze^ = I viT^^ ^^-^'^ 

is finite. At small z, one can set V{T) ~ 1, r ~ 1/^, so the contribution from this region 
is finite (when d > 2). Thus, the fate of the massless meson lies in the large z region. We 
will return to it in later sections and see that for sensible potentials, the normalization 
integral (|3.25 ) is divergent, so there is no massless scalar in this sector. 



So far, we have discussed the physics associated with the TDBI action of a real tachyon 
) . In QCD the tachyon is complex, and it is natural to study the action (|3.2|) describing 



a complex tachyon and gauge fields A^^\ A^^^ that couple to the chiral currents. In this 
case it is natural to write the complex tachyon field as T = rexp(z^), and define vector 
and axial gauge fields V = A^^^ + A^^\ A = A^^^ — A^^\ The vacuum solution can be 
taken to he 9 = 0, V = A = 0, and r = r(r). Then the complex action ( |3.2| ) reduces to 
the TDBI action for a real tachyon field, and r(r) coincides with To(r) defined earlier in 
this section. 

In this system, there are four types of small fiuctuations which give rise to mesons. 
Scalar and vector mesons correspond to fiuctuations of the fields r and V, and are governed 
by the same equations as in the real TDBI case. Axial vector mesons arise from transverse 
fiuctuations of the gauge field A, and are described by the quadratic action 

S2^-J d'^+'xV^V{r) (^G'^'^'c'^'^'F^^lFiX + G^^t^AmA^) • (3.26) 
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The first term is the same as in ( |3.22|) (up to a coefficient from the way we defined the 
vector and axial fields here). The second is due to the fact that T is charged under the 
axial gauge field. Following the same steps as above (which are discussed in detail in 
Appendix A), we find that axial vector mesons are obtained by solving the Schrodinger 
problem (|3.19|) with the effective potential 



yi^) = i(i?r- V + (2rr(r))2. (3.27) 



Here B{z) is given by (|3.24|) , and r{z) is the solution of (|3.16|) , as before. The difference 
between the vector and axial case lies solely in the last term in ( |3.27| ) . It is negligible in the 
UV, but in general important in the IR. Pseudoscalars arise from longitudinal fluctuations 
of the axial gauge field and fluctuations of 6*, the phase of T. We discuss them in Appendix 
A. 

In this section, we presented the general structure of the vacuum and excitations of 
the tachyon DBI system in the vicinity of the CPT. We next analyze some aspects of the 
dynamics for different classes of potentials. We start in the next section with a discussion 
of potentials V{T) which vanish at a finite value of T, and then move on to those that 
asymptote to zero as T — )■ oo. 



4. Hard wall models 

In this section, we consider models in which Tm is finite and the tachyon potential 
V{T) takes the form in figure 4(a). We will refer to such models as hard wall models, for 
reasons that will become clear shortly. We assume that as T — )■ Tjr, the potential behaves 
like 

V{T):^iTiR-Tr (4.1) 

where a is a positive real number. This includes the defect fermion system of section 2, 
but is of course not the most general behavior. It is easy to generalize our discussion to 
other cases, but we will not do so here. 

Plugging ( [4.1| ) in the equation of motion (|3.6|), one can check that any solution that 
reaches T = Tjr behaves in its vicinity like 

T7fl-To(r)^cv/F^; c=^^^^tl), (4.2) 
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where ^ is a dynamical scale that can be calculated by matching the behavior near Ti^ to 
the behavior ( |3.8| ) near the UV cutoff. Note that the solution ( [4 .21) satisfies the boundary 
condition (|2.13| ). Indeed, the function r(T) has a quadratic minimum at Tj^, r(T) — /i = 
(T — TifiY /c^. As T varies between zero and Tj^, r varies between A and ii. Thus, 
fluctuations described by the TDBI action see a space cut off in the IR, as in hard-wall 
AdS/QCD and the models of fg^fg . 



To study the spectrum of mesons, it is convenient to change coordinates from r to z 
using (|3.16| ). Since 2 is a decreasing function of r, an important question is whether the 
region r > |U is mapped to a compact region in z. It is easy to see that the answer is 
affirmative. Indeed, when r c:^ fj,, one can replace the factors of r in ( p.l6|) by n and use 
the fact that TQ(r) — )■ 00 as r — to rewrite (|3.16| ) as 

-^{^z-To) = 0. (4.3) 

Since To(r) approaches a finite value Tm as r — /x, z{r) must approach a finite value zjn 
as well. To calculate this value, we need to solve the full e.o.m. ( ^.6| ) and integrate ( |3.16| ) , 
which is typically only possible numerically. 

The discussion of the last paragraph is valid for any model with finite T/jj, fj,. Speci- 
fying to potentials of the form (ffTl) and using (^^), we find that 



ZiR - z{r) = -y/r - 11. (4.4) 

In section 3 we saw that scalar perturbations are described by the Klein-Gordon Lagrangian 
( ^.13| ) with a non-trivial metric ( |3.12[ ), dilaton and mass function ( ^.14| ). For potentials 
V{T) that behave like ( [4.1[ ) near T//^, we can calculate the form of the dilaton ^{z) and 
mass function m{z) near zjr. They turn out to be 

$(z) =-{a + 2)\iv{ziR- z), 

2 a (4.5) 

IJ-^ziR - zy. 

Plugging eqs. ([4.4|), (|4.5|) into (|3.20| ), we find that the effective potential for cr- mesons is 



Ves{z) - ^- '- ^ 4.6 

4 [ziR - zy 

near z = zjr. The spectrum of cr-mesons is obtained by solving the Schrodinger equation 
( p.l9|) , with a potential that behaves like (|3.21| ) for small z and ( [4.6| ) for z ~ zjr. The 
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boundary conditions for the wavefunction t/j are Dirichlet at 2; = and Neumann at 

Z = ZiR. 

It is interesting that the potential Veff diverges like A/ [z — near both boundaries 
Zb = 0, z/fj. As is well known (see e.g. [0), for such potentials the coefficient A is 
important. Near z = 0, this coefficient takes the critical value A = —1/4 separating 
different regimes; near zm it is at the critical value when a = 1 and otherwise strictly 
above it. In fact, for a > 2 the potential goes to infinity as ^ — j- zir. 

To solve for the full spectrum of a-mesons requires knowledge of the full potential 14ff 
which interpolates between the asymptotic behaviors ( |3.2ip and ( [1.6| ). However, for highly 
excited states one can neglect the potential, and read off the spectrum from the locations 
of the walls. This gives 

m„ ~ . (4.7) 

ZiR 



For the defect example discussed in and section 2, one can solve (|3.16| ) numerically, and 



show that ziR ~ 1.6//U. Plugging into ( [4.7|) gives the asymptotic spectrum 

rrin ~ 1.96/in, (4.8) 

in good agreement with eq. (4.11) in 0, m„ ~ 1.97/un. 

It is easy to repeat this discussion for vector mesons. The dilaton is given in this case 
by ( |3.23[ ), which behaves near zir as 



^{z) = -a\iv{ziR-z). (4.9) 

Plugging into ( p.20|) (with m(r) = 0), we find that the effective potential for vector mesons 
behaves near zir like ( [4.6|) as well. Hence, we expect the spectrum of vector mesons to 
have the same asymptotics as the scalar spectrum. Eq. (4.22) in gives rrin — 1.96/un, 
in excellent agreement with ( [4.8| ) .B 

In section 3, we mentioned that longitudinal fluctuations of the vector field can in 
principle support a massless scalar meson if the integral A^^, ( |3.25D is finite. The only 
possible divergence of this integral comes from the region z — )■ zir. Plugging (^7T|), ([4.2|), 
and (U) into ( ^) , we find 

/""^^ dz 



''' The agreement for vectors is slightly better than that for scalars since in Q we deduced the 
asymptotic spectrum of vector mesons from more states than for scalars. 
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It is finite for a < 1, and divergent otherwise. This is a reasonable result, since for a < 1 
V'{T) diverges as T — t- Tm. If we demand that the gradient of V remains finite, the 
massless scalar excitation is non-normalizable. 

To summarize, we find that any system described by the TDBI action ( |2.8|) , with 
a tachyon potential V{T) that behaves as with m slightly below the BF bound 

for small T, and like ( [4 .11 ) near the point Tm where the potential vanishes, dynamically 
breaks conformal symmetry. The vacuum is described by a configuration T = T'o(r), which 
behaves like ( |3.8|) , ( |3.1C1| ) for large r and like ( f4.2|) near Tjn. The spectrum of scalar and 
vector mesons is discrete and has the asymptotic Kaluza-Klein form ( [4.7| ). Although these 
results were obtained for potentials with the asymptotic behavior ( [4.1| ), we believe that 
the qualitative picture holds much more generally for models with finite T/^. We will leave 
further study of such models to future work. 



5. Soft wall models 

The models described in the previous section have the property that the asymptotic 
high mass spectrum has the KK form ( [4.7| ). In some applications, this is an undesirable 
feature. In particular, in large QCD it is believed that the masses of highly excited 
mesons behave asymptotically like 

ml ~ n, (5.1) 

sometimes referred to as a linear confining spectrum. Hence, it is of interest to try to 
accommodate such behavior in a holographic framework.i In AdS/QCD, this was achieved 
by replacing the hard IR cutoff on the radial coordinate with a soft cutoff corresponding to 
a non-trivial radial dilaton profile, chosen to get linear confinement More elaborate 
constructions appear in [^ ,| 



In this section, we will see that linear confinement and other non-trivial behaviors 
can be achieved in our framework by considering potentials V{T) with Tjji — )■ oo. Such 
potentials have the qualitative form of figure 4(b), and appear naturally in the study of 
open string tachyon condensation in string theory [|T6| - |T9| , |23| - |25[| . We will refer to them as 
soft wall potentials. 



^ Of course, near the CPT of QCD there is no particular reason to expect the highly excited 
meson spectrum to have the form ( ^.l]) , but it may play a role in this and other applications. 
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As in the previous section, we will not analyze the most general possibility. As a first 
step, we consider the class of potentials that behave at large T like 



y(T)~e-^^'^' (5.2) 

where /? is a positive constant. Potentials of this form (with a fixed value of (3) naturally 
arise in studies of open string tachyon condensation ||4^,0 . There is no particular reason 
to expect them to play a role in AdS/QCD, but we will see that they give rise to interesting, 
/3-dependent, physics. 

Proceeding as in section 4, we substitute the potential (|5.2D into the equation of motion 
( p.6|) to find the form of the vacuum solution To(r). The small r behavior of this solution 

To(r)^f^)\ (5.3) 



, r . 

It has the qualitative form depicted in figure 5, in agreement with the discussion of section 
3. In particular, there is no excised region in r, like in soft wall AdS/QCD pi]] . 

As before, we change variables from r to z using the map ( |3.16| ). At small r, this map 
reduces to ^ 

<r)--{^Y. (5.4) 



Note that as r — )■ 0, ^ — )■ oo, so the z coordinate runs from zero to infinity as well. 

We can again compute the effective dilaton and mass function for scalar mesons ( |3.14| ), 
which (for large z) take the form 

$(2),m2(z) ~ ~ (^z)^. (5.5) 

Plugging into ( |3.2C1| ), we find the large z behavior of the effective potential. 



Ves ^'{^zy—^ . (5.6) 



When /? < (i, the effective potential (|5.6|) goes to zero as z — )■ oo. Hence the spectrum 
contains a continuum of "scattering" states living at large z. When (3 > d, the potential 
grows without bound as ^ — )■ oo, and one expects a discrete spectrum. The asymptotic 
form of the spectrum is determined by the behavior of the potential at large z. One can 
read off the asymptotic spectrum from a similar analysis performed in : 



Mi-I). (5.7) 
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The above discussion was for scalar mesons, but it is easy to check that the behavior 
( ^■6|) , (|5.7|) is similar for vector mesons. The two cases differ only in the proportionality 
constant in ( |5.7| ) . This difference comes from the last term in the effective potential ( |3.20| ) 
for scalars, which contributes to the coefficient and is absent for vectors. 

In section 3 we mentioned the possibility that longitudinal vector excitations give rise 
to a massless scalar meson, depending on the finiteness of the normalization integral ( |3.25|) . 
It is easy to see that for potentials with the behavior ( p.2| ) this integral is always divergent, 
so such a scalar is absent. 

For axial vectors, the effective potential is given by ( |3.27|) , which differs from the 
effective potential for vectors by the last term ~ (rr)^. Using ( ^^) , ( |5.5|) it is easy to 
check that this term contributes to the coefficient of 22(/3-d)/(/3+d) ( |5.6[ ). Thus, the 
asymptotic spectra of vector and axial vector mesons in models with potentials of the 
form (|5.2| ) are different.! This is at first sight puzzling, since one might think that large 
energies correspond to large r, or small z, where the background r(r) goes to zero and 
chiral symmetry of the Lagrangian ( p.2| ) appears to be restored. We will return to this 
issue and explain why this intuition fails in section 7. 

We see that soft wall tachyon DBI models naturally give spectra different from those 
of hard wall models. Note that as /3 — oo, i f).7\) approaches a linearly confining spectrum. 
Thus, it is natural to ask whether there exists a potential that reproduces the behavior 
( pTT|) exactly. Presumably, as T — )■ cxd such a potential would have to go to zero faster than 
( ^.2|) with any finite (3. 

With this in mind, we next consider a potential V{T) with the large T asymptotic 
behavior 



Here 7 is a positive constant whose value will play a role below. As is by now familiar, we 
start by solving eq. ( p.6|) for the large T form of the vacuum configuration T'o(r), which 
turns out to be 



V{T) ~ e 



—jT'^ InT 





(5.9) 



We make the change of variable ( |3.16| ) and find that in the IR 




(5.10) 



The same is true for scalars and pseudoscalars. 
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The inverse map is r ^ |u(ln ^z) ^Z^'''. The dilaton behaves at large z hke $ ^ — In V{Tq) = 
-fT^lnTo. Using eq. OCT}] ) gives 

^{z) c^-fH^z'^{lnz)^~^. (5.11) 

When 7 = (i, one finds $ ~ /U^-z^; the potential KfJ ( |3.2C1| ) for scalar mesons and its analog 
for vector mesons goes in this case like /U^-z^, which leads to the linearly confining spectrum 
~ n'^n. For other values of 7 the quadratic potential receives logarithmic corrections, 
which deform the linear spectrum (|5.1D . 

The coefficient of fi'^n in is different for scalar and vector mesons. This is because 
the mass term in the effective potential (the last term in (|3.20| )) contributes to the coeffi- 
cient of z'^ in Veff for scalar mesons, and is absent for vector mesons. On the other hand, 
the difference between the highly excited spectra of vector and axial vector mesons lies the 
last term of the axial vector potential ( p.27|) , which goes like (rr)^ ~ (r^^)^ ~ z^/(lnz)^. 



It is (marginally) subleading at large z, so the spectra of vectors and axial vectors (slowly) 
approach each other at high excitation levels. 

To summarize, in this section we saw that potentials of the general form of figure 4(b) 
give qualitatively different physics than those of figure 4(a). Rather than leading to an 
excision of a finite region of small r, they provide a soft wall that smoothly suppresses 
dynamics at small r. By varying the potential, one can arrange the spectrum to have 
different asymptotic forms such as (|5.1| ) and ( p.7| ). All the features described above arise 
dynamically once a potential V{T) is specified. This is in contrast to the usual approach 
where the profile of the dilaton and other closed string fields are put in by hand. As we 
will discuss in section 7, in QCD the difference between the two approaches is that one is 
obtained by expanding around the boundary between confining and conformal behavior, 
while the other is obtained by expanding around the theory with very few flavors. 



6. Finite temperature 

In this section we will study the finite temperature thermodynamics of systems de- 
scribed by the TDBI action (pT^). We introduce temperature in the standard way, by re- 
placing the background AdSd+i in which the tachyon field lives with an AdS/Schwarzschild 
black hole which acts as a heat bath. The metric is given by 

ds' = (^-Fir)dt' + B^^O'] + ^r) = 1 " (y)' (6-1) 
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where = tcT is the location of the horizon. T here is the temperature, not to be confused 
with the tachyon field. To avoid any confusion, in the rest of this section we will use rh to 
label the temperature. 

We note in passing that in ^ the temperature was introduced in a slightly different 
way. There, the background in which the tachyon lived was always AdS^ x S^, and at 
finite temperature we replaced it by an AdS^ black hole, which amounts to taking F{r) = 
1 ~ {i^h/i^)'^ for all d in (|6.1J ). Presumably, the different ways of introducing temperature 
should give similar results for the thermodynamics. 

The TDBI action (density) at finite temperature is given by 

S = - j drr'^-^V{T)^l + F{r){rdrT)^, (6.2) 
and the equation of motion which follows from it is 

r^^-i^T In y(T) d ( r'^+^FT' 



1 + Fr2T'(r)2 dr \ + Fr^T'{r)^ 



(6.3) 



As a check, for zero temperature {rh = 0), (|6.3|) reduces to (|3.5|) . It is useful to change 
variables from r to 



, d 

r 



y-[-), (6.4, 

in terms of which the action (|6.2|) takes the form 

S = ~ri j dyV{T)y/l + F{y){dyT')\ (6.5) 

where F{y) = 1 — ^ and T' = dyT. Note that in (|6.5D, the explicit dependence on rh has 
been factored out. The corresponding equation of motion. 



+ F(y){dyVf dy \ + F(v)(dyT) 

is independent of rh when written in terms of the variable y. y ranges between 1 and 
(A/r^)*^. We will mostly work in the BKT limit discussed in and the previous sections, 
in which 1 < y < oo. 

Some aspects of the finite temperature analysis are different for the hard wall (in 
which the potential V{T) vanishes at a finite value of T = Tm) and soft wall (infinite Tifi) 
models. We will start by discussing the situation for the soft wall case, and then mention 
the details that are different in the hard wall one. 
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We saw in section 3 that at zero temperature, the ground state of soft waU models 
corresponds to a tachyon configuration of the sort plotted in figure 5. To find the config- 
uration at finite temperature, we can proceed as follows. Denote by Th the value of the 
tachyon field at the horizon, 

T{y = 1) = Th. (6.7) 



Assuming that T{y) is smooth near y = 1, the equation of motion ( |6.6| ) then determines 
T'{y = 1) to be 

T'{y = l) = ^dT\nV{TH). (6.8) 

Thus, the entire solution of (|6.6| ) is determined by specifying Th. 

The value of Th is determined by the temperature rh in the following way. As y — oo 
the tachyon field goes to zero, and we can replace the full equation of motion ( |6.6D by its 
linearized version, 

^ W - ym -It = o. (6.9) 

The general solution of this equation behaves at large y as 

T{y) ~ ^(cilny + C2) (6.10) 

where the constants Ci are functions of Th. 

On the other hand, as discussed in 0, the UV boundary conditions imply that the 
asymptotic large y behavior of the finite temperature profile should be the same as at 
zero temperature. The latter is given by eq. ( |3.10|) , which can be written in terms of the 
coordinate y (up to an unimportant overall constant) as 

Toiy) ^ Iny + C2 + Ci In (6.11) 

The ratio of the coefficients of the constant and log terms in ( |6.1C1| ) and ( |6.11| ) must be the 
same 0, 

= ^ + (6.12) 

dci(Th) Ci n 

Given a potential V{T) and a definition of the scale fi, ( p.l2| ) gives a relation between 
and Th. 

A useful observation for the discussion below is that the value of the tachyon at the 
horizon, Th, goes to zero at a finite value of the temperature, rh = r^h^^^^ . To determine 
r^h"^\ one can proceed as follows. For small Th, the full solution of (|6.6| ) is in the small field 
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regime, and one can determine T{y) by solving ( |6.9| ). This is the hypergeometric equation 
with a = b = 1/2, c = 1. Its generic solution diverges (logarithmically) as y — )■ 1 and is 
thus not suitable for our purposes. The solution that satisfies the boundary conditions 
T{y = l)=n is 

T{y) = -TuK{l - y) (6.13) 

TT 

where 2K{x)/ti is the hypergeometric series 



K{x) is known as the complete elliptic integral of the first kind. As y — )■ oo, it behaves like 

i^(l - y) ^ Q Iny + 21n2 + ■ ■ (6.15) 
Comparing to ( |6.10| ), we see that as Th — )■ 0, the ratio C2/C1 approaches the finite limit 

lim =41n2. (6.16) 

TH^QCiiTh) 

Plugging into ( |6.12| ), we find the critical temperature 



rr^^=2^/.exp(^-^J, (6.17) 

which is given in terms of quantities characterizing the zero temperature problem, namely 
H and Ci. 

As Th approaches r]^^^^\ goes to zero. In this limit, the free energy of the non-trivial 



solution of ( |S.3| ) approaches that of the trivial solution T = 0. In general, the difference 
between the two is given (up to an overall multiplicative constant) by 

/•A 

J^{rh)= / drr'^-^ V{T)^/l + F{r){rdrT)^ - 1 . (6.18) 

Here we have returned to the r parametrization, and are working at a large but finite cutoff 
A, in order not to have to worry about UV divergences. The BKT limit can be taken at 
the end of the calculation. 

We would like to evaluate the derivative of J-' w.r.t. r^. For this purpose we consider 
the quantity 

SJ^ = J^{rh + SrH)-T{rH). (6.19) 
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This difference receives contributions from two sources: 

(1) Tlie range of tlie r integral is different in tlie two terms. Tlius, one contribution to 
tlie difference is 



drr 



d-i 



V{T) + F{r){rdrT)^ - 1 



(6.20) 



To first order in 5rh we can evaluate the integrand of ( |6.2C1| ) at r = rh, and multiply 
hj Srh, yielding 

dJ'^'') = drHri-\l-V{n)) (6.21) 



where we used the fact that F(rh) = and drT{r = rh) is finite (see eq. (|6.8|)). 
(2) When we change rh by 5r^, the whole solution T(r) changes to T{r) + 5T{r). If we 
denote the integrand in ( |6.18| ) by C{T,T'), the second contribution to ( |6.19| ) is 

dr [C{T + 5T, T' + 6T') - C{T, T')] . (6.22) 

rh+Srh 



Since 5T is of order 5rh, we can replace the lower limit of integration in ( |6.2C1| ) by rh, 
and write 



dr 



5C 5C ■ 



5C 



ST 



A 



(6.23) 



In the last step, we integrated the second term by parts and used the equation of 
motion of T. The contribution from r = A vanishes since dT is zero there. The one 



from r = Th vanishes as well, since is proportional to F{rh), which vanishes. Thus, 
To summarize, we conclude that 



drh 



ri-\l-V{Th)). 



(6.24) 



This is a non-trivial equation, since Th itself is a function of rh, which we have to find first, 
given the potential V{T). Once we know it, we can solve ( |6.24| ) to find J^{rh)- A general 
property we can read off from ( |6.24|) is that J-" is a monotonically increasing function of rh- 
In appendix B we numerically study the finite temperature behavior for some particular 
soft wall potentials; our results are compatible with ( |6.24|) . 



Expanding ( |6.24| ) around 



,(crit) 



( |6.17| ), we find to leading order: 



dF d^ 



drh 8 



rt'Ti + 0{T^) 



(6.25) 
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We can also compute T for small T^, directly, by expanding ( |6.18|) in a power series in T. 
To determine the form of T up to order T^, one needs to keep contributions to ( |6.18|) up 
to order T". The leading contribution is quadratic, 



(6.26) 



Integrating the first term by parts we get the e.o.m. (|6.9|) , so the integral ( |6.26| ) reduces to 
evaluating r'^+'^ F {r)TT' at the boundaries r = r^, A. Both of these vanish, since F{rh) — 
and T(A) = 0. 

Thus, the leading contribution to J-" is quartic in T^,. To calculate it we need to include 
the quartic contribution to the potential V{T), 



Vm = 1 + -m^T^ + -T^ + 



(6.27) 



Plugging ( |6.27|) into (|6.18|) , we find 



^4 



drr 



d-i 



F{r){rT'f + m^T^ + -T^ + -m^F{rTT'Y - ^F'^{rT 



/\4 



(6.28) 



The solution to the full e.o.m. ( |6.3| ) can be expanded as T{r) = Ti(r) + T3(r) + . . . where 
Ti{r) is of order T^. In particular, Ti is given by (|6.13| ), while T^ir) is obtained by varying 
( |6.28|) and keeping only terms cubic in T^. 

Using the e.o.m. from ( |6.28|) and plugging back into J-4, one finds (in the BKT limit) 



^4 



'rh 2 



(6.29) 



Since all terms in ( 6.29 ) are quartic in T, to leading order in Th we can replace T by Ti 
( |6.13| ). Upon performing the integral, we arrive at 



where 



ci ~ -2.847, C2 ~ 0.072. 
Thus, we see that the leading behavior of J-" is 



(6.30) 



(6.31) 



T = cTt + 0(T,6) 
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(6.32) 



where 

d 



c= ^(rf^^Ylcia + rf^). (6.33) 



To compare to ( |6.25|) , we need to differentiate (|6.32| ) w.r.t. r^. This gives 

AcTlT'^ + --- (6.34) 



where T'^ is the derivative w.r.t. r^. Comparing ( |6.25| ) and ( |6.34D and substituting = 
r)^ , which is okay to leading order in T^^, we find 

n = ^yr'Y-\r, - rt'^) + . . . (6.35) 

We are now ready to discuss the finite temperature phase structure of soft waU TDBI 
models. The phase structure is different for c < and c > 0, and we will discuss the two 
cases separately. 

Consider first the case c < 0. There are two regions in the phase diagram in which we 
understand the properties of the non-trivial solution T(r): low temperature [i.e. small r^), 
and the vicinity of the critical temperature ( |6.17| ). As rh — )■ 0, the value of the tachyon 
field at the horizon — )■ oo, and the free energy J-" approaches a finite negative constant. 
Near r^i^'^^^\ is small. Note that eq. ( |6.35|) implies that the non-trivial solution only 
exists for rh < t^^"^^^^ • '^^^ ^^^^ energy goes to zero as well, 

^~-(r,-rf^*^)^ (6.36) 



as can be seen by combining ( |6.32| ) and (|6.35|) . 



The most natural way to connect the two behaviors is depicted in the left two panels 
in figure 6. These figures imply that the system exhibits a second order phase transition 
at Th = ^i'^'^^*^- In the low temperature phase, the solution with lowest free energy has 
a non-trivial T{r). As rh approaches the critical value, it smoothly goes to the trivial 
solution T = 0, and above the transition temperature, only the trivial solution exists. 
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c<0 



c>0 



Fig. 6: Qualitative behavior of rh{Th) and J-irh) for c < (left) and c > (right). 

For c > the situation is the foUowing. At low temperatures (smaU rh) things are as 
before, but near = r^/^^^^^ they are different. In particular, (|6.35| ) implies in this case that 
there is a non-trivial solution with arbitrary small Th, but it exists only for rh > r^h^^^"^ ■ 
The natural way to connect the two behaviors is depicted in the right panels in figure 
6. The upper panel shows rh{Th). Starting from Th = 0, must initially increase, as 
indicated by ( |6.35| ), but eventually it must turn around and asymptote to zero as T/^ — )• oo. 
Assuming that there is a single turning point leads to the upper right plot in figure 6. 

The corresponding behavior of J^{rh) is depicted in the lower right panel of figure 6. 
At low temperature, J-" is negative and grows with rh- At a maximal value of rh, which 
can be read off the upper right panel, it connects to another branch of solutions that joins 
smoothly to the behavior near Th = ( |6.35|) . In this case, the transition is first order; 
it happens when the free energy of the low temperature branch of non-trivial solutions 
approaches that of the trivial, T = one. It occurs at the temperature at which the lower 
curve in the lower right panel in figure 6 intersects the rh axis. 

The strength of the first order phase transition is controlled by the size of the coefficient 
c ( |6.33|) , which in turn is determined by the coefficient a of the quartic term in the potential 
( |6.27|) . As c decreases, the transition becomes less strongly first order, until eventually c 
flips sign and the transition becomes second order, as discussed above. If one tunes a such 
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that c vanishes, the leading behavior of T near = is as a higher power of than 
( |6.32|) , and one must redo the analysis taking into account higher order terms in 
We will leave this analysis to future work. 

The plots of figure 6 were obtained by assuming the simplest possible behavior given 
the analytic results obtained above. To see that this is indeed what happens, we analyze in 
appendix B some particular potentials, and verify numerically that the behavior of figure 
6 is indeed reproduced, both for positive and negative c. 

Our discussion so far has concerned soft wall models. Some aspects of the analysis 
are the same in the hard wall case. The differential equation ( |6.24 ) that encodes the 



dependence of the free energy on temperature, and in particular the monotonicity of this 
function is still valid, as is the analysis of the small Th region, which only relies on the 
form of the potential ( |6.27| ) for small T . 

On the other hand, the small r/j (low temperature) analysis is different for soft and 
hard wall models. For the soft wall case, the solution T(r) ends on the horizon for any 
finite temperature, and the maximal value of T(r), which we denoted by above, is finite 
(and thus strictly below the infinite Tjji). For hard wall models, if the temperature is low 
enough the vacuum solution To(r) does not reach the horizon, and the maximal value of 
To(r) is equal to Tjji. 

The phase diagram is still expected to be different for c < and c > 0. In the former 
case, it is expected to be given by the left panel in figure 7. For small rh, the non-trivial 
solution T(r) does not intersect the horizon; its free energy is given by the blue curve. At 
some value of the temperature, corresponding to the point where the blue and red curves 
meet, the solution first touches the horizon, and above this temperature it behaves like in 
the soft wall discussion above. In particular, the value of the tachyon on the horizon is 
smaller than Tjji, and goes to zero at a critical temperature r^i^^^^\ which corresponds to 
a second order phase transition. 



I'h 'h / ""h 



c<0 c>0 

Fig. 7: Qualitative behavior of rh{Th) for hard wall models with c < (left) and 
c > (right). The red (blue) curve corresponds to solutions that do (do not) 
intersect the horizon. 
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For c > we expect a phase diagram of the sort exhibited in the right panel of figure 
7. The blue curve corresponds again to solutions that do not intersect the horizon, while 
the red curve labels solutions that do. The phase transition is in this case first order; it 
occurs at the temperature at which the blue line intersects the rh axis. The situation here 
is very similar to that analyzed in section 5 of (see figure 7 in that paper) , and we refer 
the reader there for further discussion. 

In summary, for both classes of tachyon potentials, which give rise to hard and soft 
wall models, changing the quartic term in the potential leads to an interpolation between 
first and the second order phase transitions. This is in contrast with other holographic 
models in which the transition is first order, as in models with geometric confinement ||26|| , 
probe D-branes |l5T| - |5^ and hard and soft wall models [|5^ . Second order phase transitions 
in holographic systems have been observed before in holographic superconductors . Our 
setup is different; it will be interesting to explore possible phenomenological implications. 



7. Discussion 

In this paper we studied a class of continuous phase transitions in d dimensional 
(large A^) QFT known as conformal phase transitions . We argued that in studying these 
transitions one can focus on the dynamics of the order parameter, an operator whose 
dimension A approaches d/2 near the transition. Since these transitions often occur at 
strong coupling, it is natural to use holography to study them. In the bulk description, 
this corresponds to analyzing the dynamics of a scalar field T in AdSd+i, whose mass is 
close to the BP bound. When the squared mass of T goes below the bound, the theory 
develops dynamically a small mass gap fj,. Interestingly, the universal physics near the 
transition^ is sensitive to the full non-linear Lagrangian of T. We took this Lagrangian 
to have the tachyon DBI form ( p.8| ), (|3.2| ), and analyzed the dynamics near the transition 
for different potentials V. This can be viewed as a bottom-up description, in the spirit of 
AdS/QCD and AdS/CMT studies. 

We found that the dynamics of this class of models is similar to that of hard and soft 
wall models of AdS/QCD pi]|j2I[] , with a dynamically generated wall. Hard wall models 
are obtained when the potential V{T) vanishes at a finite value of T, Tjji, while soft wall 

I.e. the physics at the scale fi, which goes to zero at the transition and hence can be taken 
parametrically smaller than all other scales in the system. 
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ones correspond to potentials with Tjji — t- oo. For hard waU models, the spectrum of 
small excitations (mesons) near the transition was found to take the standard Kaluza- 
Klein form (|4.7D at large excitation level, while for soft wall ones we exhibited potentials 



with non-trivial high energy spectra (|5.1|) , (^77|). We also studied the finite temperature 
thermodynamics of these models and found that they exhibit first or second order phase 
transitions, depending on the potential V{T). Interestingly, the high energy spectrum of 
mesons and the type of finite temperature phase transition exhibited by the model are 
sensitive to different features of the potential. The spectrum is sensitive to the behavior 
of V{T) near the bottom (see figure 4), while the finite temperature phase transition is 
sensitive to the behavior near the maximum at the origin of field space. 

In this section we would like to comment on some aspects of our construction and 
possible generalizations. One natural question is why it is sensible to focus on the dynamics 
of the bulk tachyon field T dual to the boundary operator O whose dimension approaches 
d/2 near the transition, and neglect the backreaction of other fields on it. The reason is 
that near the transition one can think of the boundary theory as a CFT deformed by the 
"double trace" operator O"^. Correlation functions in the deformed theory, such as the two 
point function (((9(x)(9(0))), are given by expressions like 



((C(a;)C(0))) = (C(x)C(0)exp 



-A j d'^zO^iz) 



) (7.1) 



where the expectation value on the right hand side is computed in the CFT. Thus, the 
two point function ( [7.1|) can be calculated in conformal perturbation theory in terms of 
(integrated) correlation functions involving multiple O operators in the CFT. 

The dual statement in the bulk theory is that in order to calculate properties of the 
massive theory, all we need to know is the non-linear Lagrangian for the order parameter 
T{x^). To compute correlation functions of other operators in the boundary theory, such 
as the global symmetry currents discussed above, we similarly need to know the non-linear 
terms in the Lagrangian involving the corresponding bulk fields and an arbitrary number 
of T's, which is precisely the type of information encoded in actions such as ( |3.2| ). Of 
course, the above discussion does not imply that the dynamics of T should necessarily 
be governed by the TDBI action with any potential. The assumption that it is, is the 
bottom-up aspect of our analysis. It would be nice to provide further justification for this 
assumption or to improve on it. 

It is important to stress that our construction is quite different from what is known 
in the study of holographic systems as the probe approximation. In the context of QCD, 
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this approximation is valid in the 't Hooft hmit, in which the number of colors goes to 
infinity while the number of fiavors F is kept fixed. In this limit, the leading contribution 
to the bulk action involves the metric and other closed string fields, and takes the form of 
Einstein gravity with various corrections. Confinement should be a property of this action, 
and manifest itself via the generation of a hard or soft wall in the geometry. 

The first subleading term in F/N involves fields associated with the quarks, such 
as the tachyon which is dual to the meson operator, the gauge fields dual to the global 
currents etc. Their action is usually taken to be of the TDBI form ( |3.2|) with corrections 
(see e.g. [^^,[|5^). Although this term is subleading in F/N, it is the leading contribution 
to the dynamics of fields in the meson sector, and clearly needs to be kept when studying 
meson dynamics. 

As F/N increases, the probe approximation becomes worse and worse, until eventually 
when F ~ A^, it breaks down. In general, in this region of parameter space we have no 
reliable tools for studying the dynamics for energies of order Aqcd and below, where the 
theory is strongly coupled. The main point of our construction is that as F — )■ Fc {i.e. 
near the transition from conformal to confining behavior), we again expect to be able to 
provide a good bulk description of the dynamics. However, the small parameter is now 
not F/N, as in the probe approximation, but the ratio of scales ^/ Kqcd-, where ^ is the 
dynamically generated meson mass scale, and Aqcd is the crossover scale between the free 
UV and interacting IR CFTs. 

As we argued, near the transition we can describe the dynamics by writing an action 
for the bulk field T dual to the meson operator. However, this action plays a very different 
role from that which appears as the first subleading contribution to the bulk action in the 
probe region F N. In particular, we should not add to the bulk Lagrangian an Einstein 
term for the metric and dilaton fields, which would drastically modify the dynamics of 
the tachyon. The interpretation of this modification in the boundary theory would be 
that in addition to the operator (9^, the Lagrangian includes non-trivial couplings of other 
marginal or relevant operators, which is not expected to be the case near the transition. 

Another interesting issue concerns the chiral symmetry of the complex TDBI action 
( pi.2| ). In the massive phase, chiral symmetry is broken by the non-trivial condensate of 
the tachyon field T. This condensate goes to zero at large r (or small z); therefore, the 
symmetry is restored there. In the boundary theory, the restoration is refiected in the 
short distance structure of off-shell Green functions of the currents dual to the gauge fields 



36 



Naively, one might expect the restoration to also be visible in the highly excited 
spectrum of mesons. However, we found by an explicit calculation that for potentials 
of the form ( p.2| ) (at large T), the vector and axial vector spectra remain distinct at 
arbitrarily large excitation levels-IHI This is not surprising given that the highly excited 
spectrum depends on the large z (or small r) behavior of the vacuum solution T'o(r). This 
region is usually thought of as the IR region in terms of the usual correspondence between 
scale in the field theory and position in the radial direction in AdS^ and in it the symmetry 
breaking effects due to the non-zero To(r) are not suppressed. 

At first sight it is surprising that the masses of highly excited mesons, which are 
very heavy, are sensitive to IR physics. In fact, this is a phenomenon familiar from black 
hole physics and holography, known as the UV/IR correspondence. The TDBI analysis 
incorporates the fact that highly excited mesons are physically large (in IR*^), hence their 
properties are sensitive to the long distance behavior of the theory. In particular, there is 
no contradiction between the fact that chiral symmetry is restored at short distances in off- 
shell Green functions, and its lack of restoration in the highly excited meson spectrum (see 
|5^,^ and references therein for further discussion). In fact, the TDBI system provides a 
nice laboratory for studying the relation between the two in a large class of models. 

One of the important applications of strongly coupled models of the sort we studied 
is to electroweak symmetry breaking (technicolor). In this context, the system analyzed in 
this paper can be viewed as the symmetry breaking sector of walking technicolor. We will 
leave a more detailed investigation of the resulting phenomenology to future work. One 
aspect that we would like to mention is the techni-dilaton, a pseudo Nambu-Goldstone 
boson of broken conformal symmetry, which may play an important role in current collider 
experiments (see e.g. for a recent discussion). 

In our model, the dilaton is the lightest cr-meson, i.e. the lightest normalizable eigen- 
state of the bound state potential ( p.l9| ), ( p.20| ). A natural question is whether this state 
is expected to be anomalously light, compared to other scalar and vector mesons. We have 
not studied this problem in our paper, but qualitatively expect this to be the case, for the 
following reason. 

As discussed in section 3, the effective potential for cr-mesons ( |3.20| ) behaves at small 
z like ( |3.21| ), and as explained there, if this was the exact behavior of the potential, the 



A similar result was found in [48|. However, we also found an example ( ^.81) where the highly 
excited spectra agree, which was not the case there. 
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spectrum would contain tachyons with mass squared of order —fi'^. The potential ( 3.20| ) 
actually deviates from ( |3.21| ) at 2; ~ l/^u and larger, which raises the mass squared of the 
lowest lying a-meson by an amount proportional to n^. Thus, it is natural to expect the 
mass of the lowest lying meson to be of order Ai^ = a^^ with a smaller than the typical 
separation between states. In we analyzed a particular (hard wall) model and found 
that this is indeed the case. 

Note that a cannot be negative, since that would imply that the configuration T = 
To(r) is unstable, contrary to assumption. For generic potential V{T) we would expect it 
to be non-zero. We do not know whether it is possible to tune the potential such that a 
vanishes, giving a massless techni-dilaton. Since the breaking of conformal symmetry in 
our model is explicit, we would expect this not to be the case. It would be interesting to 
see how small one can make a by tuning the potential. 

Another natural question is whether one can study meson dynamics in QCD using our 
approach. Our model is rather different from other holographic QCD models, primarily 
since the latter attempts to describe the physics in the limit N ^ 00, F fixed and then 
continue to the physical regime, while our approach expands around the boundary between 
confining and conformal behavior F = Fc- The two expansions can potentially provide 
complementary information about the physics of mesons. 

In Section 6, we studied the dynamics of our model at finite temperature. We found 
that the phase diagram of the system can exhibit either first or second order phase transi- 
tions, depending on the value of the quartic term in the tachyon potential. Recently there 
has been a lot of discussion of possible applications of holography to condensed matter 



physics (see e.g. ||6T| - |65[] for reviews). Second order phase transitions might be interesting 
in this context, since near such transitions one may hope to extract universal quantities, 
such as critical exponents. For example, from ( |6.36| ) one can infer the behavior of the heat 
capacity near the phase transition: it experiences a jump from a constant non-zero value 
to zero, which is the standard mean field behavior. One can tune the tachyon potential 
such that the critical exponent a associated with the specific heat obtains a non mean field 
value. We will leave a more detailed discussion of this to future work. 

Another interesting feature of our model is its low temperature behavior. In the usual 
hard and soft wall AdS/QCD models, the low temperature phase is obtained from the 
zero temperature geometry by compactifying Euclidean time. The black hole phase is 



thermodynamically favored only above a critical temperatures (see e.g. [^5|). This implies 
that there is no dissipation at small temperatures, at least in the classical approximation 
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to the bulk theory. In contrast, in our soft waU models, the black hole phase describes 
the thermodynamics at arbitrarily low temperature, i.e. the finite temperature vacuum 
solution To(r) extends all the way to the horizon at r = rh. Hence, the effective met- 
ric for fluctuations has a horizon, and the system exhibits dissipation at arbitrarily low 
temperature. A detailed investigation of this will also be left for future work. 
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Appendix A. Meson spectrum in the axial sector 

In this appendix, we discuss the meson spectrum in the axial sector. Recall that the 
action for the complex tachyon DBI ( p.2|) is 

5 = - y d'^+^xV{T^T) (y-G(^) + aZ-GW) (A.l) 

where 

^MN = 9mn + -D(^mT^ Dn)T + F^j^Ij (A.2) 

and similarly for L ■(r^ R. The complex tachyon, vector, and axial vector gauge fields are 
defined to be 

T = rexp(z^), V = A^^^ + A^^\ A = A^^^-A^^\ (A.3) 
Expanding to quadratic order, the action for the axial sector is 

S = -j d''+'xVGV{r) QG^^'C^^'Filiri^l^iv' + G'^^r2(aM^ + AM^d^e + A^)^ . 

(A.4) 

It is convenient to change to the z coordinate following (|3.16| ) , in terms of which the action 
reads 

(A.5) 
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where indices are raised and lowered with the flat metric rjap- This Lagrangian gives rise 
to axial vector and pseudoscalar mesons; we will now analyze them in turn, starting with 
the vectors. 

For the purpose of this discussion we can set 9 = 0, pick the gauge = and 
parametrize the remaining components as Aa{z,x) = Aa{x)v{z). The transversality con- 
dition takes the form = 0. The equation of motion following from the action ([A .51 ) 
is 

1 „ 1 



^d,{V{T)r{zf-^d,v{z)) + -mlV{T)r{z)'^-\{z) - V{T)r{zy-^T^v{z) = 0. 
We make the change of variable iIj{z) = e~^^'^v{z) where 

B{z) = -lny(r) -{d-3) \nr{z). 
Then ( [A .61) takes the Schrodinger form with 

V,s = \{Br-lB" + 4r{zfr{r{z)f. 
Plugging Aa{z, x) into the bulk action ( |A.4|) , we find a tower of axial vectors 

S = -jd''x + mlA^A'^^ 

subject to the following normalization conditions: 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



/ 



dzV{T)r{zY-^v{zf 
dzV{T) (r(z)^-3(t;')2 + Ar{zY-^T''v{zf) 



dzt(j{z) ~ 1 



dz 



B' 



(A.IO) 

where we have omitted unimportant overall constants. In the second line, we integrated 
by parts and used the Schrodinger equation. Since the norm is the standard Schrodinger 
one, the spectrum is obtained by solving the bound state problem for a particle moving in 
the potential ( |A.8|) . 

Next, we turn to the pseudoscalar modes, and show that the spectrum always contains 
a massless pion. Parametrizing Aa{z,x) = Lp{z)daTT{x),9{z,x) = 9{z)tt{x), the equations 
of motion of (|A.5| ) are 



^d.iViTyizf-^'d.ifi) - V{T)r{zy-^T\9 + V?) = 

T^d,9 +]m'^r{z)-^d,^ = 0. 



(A.ll) 
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The normalization conditions for the kinetic and mass terms of the field tt are 

dz Qy(r)r(2)^-3((^')2 + V{T)r{z)'^-^T^{e + ipf^ = 1 

' dzV{T)r{z)'^-^T^{djf = ml. 



(A.12) 



When = 0, the second equation in (|A.11|) implies that 9{z) is a constant, which can be 



set to zero by a gauge transformation. At small z (in the UV), the remaining equation of 
motion reduces to 

^a,(;23-d5,(^)-z(log^)V = 0. (A.13) 

The second term is presumably subleading to the first and we drop it. Then the two 
independent solutions are ~ 2;^~^,c. Plugging into ( |A.12|) , it is easy to check that the 
kinetic normalization condition does not diverge for either solution at small z, as long as 
d > 2. Since both UV solutions are normalizable, the massless mode exists as long as there 
is one normalizable solution in the IR. 

Let us confirm that this is true for the hard-wall and soft-wall models discussed in 
sections 4 and 5. For the hard wall model ( [4.1| ) at large 2, r ~ |U and r ~ Tm are both 
finite and have no effect on the divergence structure, so we ignore them here. The eom 
(up to constants) is 

d.iiziR - z^d^if) - [ziR - z^if = 0. (A.14) 

We can neglect the second term (which turns out to be self-consistent) . Then the solutions 
are ip ~ {zir—zY~'^, c. Since we are concerned with a > 1, there indeed is one normalizable 
and one non-normalizable solution. 

For the soft wall case, we first consider the potential (pT2|), V{T) ~ e~^^^ . We 
assume that {3 > d because the alternative gives rise to a continuum of scattering states 
and is uninteresting for our purposes. The equation of motion in the IR is 

(1 2^ \ 2(g-d) / 1 2/3 \ 

exp(--/3z^)a^<^j = z^+^exp ( j (A.15) 

1 2/3 

With the ansatz (p{z) = f{z)exp{-^l3zf^+''-), the eom reduces to 

f" + J^z^f' = z^f (A.16) 
p + d 
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(where again, we take a large z limit at every step). Now, we make the ansatz / 



2^ 

exp(a2;^+'* ). This yields a quadratic equation for the coefficient a: 

4a2/32 + 2a/3^ - + = 0, (A.17) 

which has two solutions, one with a > and one with a < — /3/2. One can readily check 
that the latter is normalizable. The model which gives rise to linear confinement ( p. 8] ) can 
be checked using the same steps and behaves similarly. 

Finally, we examine the spectrum of massive pseudoscalar modes. If vn? is not zero, 
we can solve the first equation in ( |A.llj ) for Q and substitute the result into the second 
one. This leads to 

where $ = V {T^r'^~'^ d z'^ . We make the additional change of variable '^{z) = e~^^'^^{z) 
where 

B{z) = InV {T)r{zy-^T{r{z)f. (A.19) 



Then the eom ( |A.18 ) takes the Schrodinger form with the effective potential 



Fefi = liB'f - + 4r{zrr{r{z)r. (A.20) 

However, the normalization condition from the kinetic term is not the usual Schrodinger 
norm. Instead, it is 



Appendix B. Thermodynamics of a particular TDBI model 

In this appendix we analyze the thermodynamics of a particular class of soft wall 
potentials in (i = 4, as a check on the general discussion of section 6. These potentials are 
given by 

The quadratic term in the potential is tuned to the BKT limit for (1 = 4, i.e.V{T) = 
1 — 2T'^ + . . . The coefficient A controls the value of the quartic term in the potential. 
Comparing to (|6.27| ), we see that a = 8 — 2 A. 
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Our goal is to analyze the thermodynamics for two values of A that give rise to 
positive and negative c in (|6.33 ). For d = 4 the transition between the two regimes occurs 
at a ^ 6.5. Thus, we will analyze two cases: 

(1) A = 0, a = 8: ( |6.33|) gives c < in this case, so we expect a second order phase 
transition. 

(2) A = 5, a = —2: c > 0, first order phase transition. 

The numerical results for A = are shown in figures 8, 9. We see that they are in 
agreement with the discussion of section 6, and in particular the two left panels in figure 
6. 

The numerical results for A = 5 are exhibited in figure 10, 11. They agree with those 
in the two right panels of figure 6, and imply that the transition is first order in this case. 



1.6 h 




Fig. 8: r^Z/i as a function of Th for A = 




-0.10 



-0.15 - 



-0.20 - 



-0.25 



Fig. 9: T I as a function of r^Z/x for ^ = 
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Fig. 10: r/i/^ as a function of Th for A = 5 



-0.005 




Fig. 11: as a function of r/i//^ for A = 5 
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